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Magnetization processes of frustrated Kondo necklace model are studied by means of a density 
matrix renormalization group (DMRG) method and an elementary band theory based on a 
bond-operator formalism. The DMRG calculations clearly show the cusp singularity in a low- 
magnetization region (0 < m < 1/2) besides that in a high-magnetization region (1/2 < m < 1) 
which is expected from previous studies on the magnetization curve of the Majumdar-Ghosh 
model. An appearance mechanism of the low-magnetization cusp is interpreted in terms of a 
double-well shape of a low-energy band arising from frustrations between nearest- and next- 
nearest-neighbor interactions. We also discuss critical behaviors of magnetization near the cusp 
and obtain a phase diagram showing whether the cusp appears in the magnetization curve or 
not. 
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1. Introduction 

Over the last few years, a large number of materials 
exhibiting quantum phase transitions induced by an ex- 
ternal magnetic field have been discovered, and their crit- 
ical properties have been attracted much attention ex- 
perimentally as well as theoretically. Such field-induced 
phase transitions are associated with various singular- 
ities in the magnetization curve, such as plateaus, 1_4 ' 
cusps 5-11 ) and jumps. 12 ~ 14 ' In contrast to concern about 
plateaus and jumps, very few attempts have been made 
at the problem of cusp singularity in the magnetization 
curve. The cusp singularity was first found in the SU(3) 
integrable spin chain known as the Lai-Sutherland model 
by Parkinson with the use of the Bethe ansatz method. 5 ' 
Subsequently, it has been investigated with a focus on in- 
tegrable models. This may be because it is difficult to find 
cusp singularities of non-integrable models by conven- 
tional numerical approaches like an exact-diagonalization 
and quantum Monte Cairo methods. 

Very recently, Okunish et al. showed that the cusp 
singularities appear in the magnetization curve of non- 
integrable spin systems, such as the spin-1/2 zig-zag lad- 
der and the spin-1 bilinear-biquadratic chain, using the 
product-wave- function renormalization group (PWFRG) 
method, 15,16 ' which is a variant of the density-matrix 
renormalization group (DMRG) method. 17 ' They also 
pointed out that the appearance of cusps originate in the 
double-well shape of the energy dispersion of the quasi- 
particles such as magnon, spinon excitations. 

In this paper, we introduce frustrated Kondo necklace 
model, which is equivalent to the spin-1/2 chain with 
nearest- and next-nearest-neighbor interactions (Ji and 
J2) coupled to a localized spin on each site, as a new non- 
integrable spin system exhibiting the cusp singularities 
in its magnetization curve. Here, let us briefly review the 
magnetization process of the frustrated Kondo necklace 
model already understood by previous studies. 18 ~ 20 ) The 
low-energy excitations of this system have a finite gap 



(so-called Kondo singlet gap), reflecting that the ground 
state in the absence of the magnetic field is always in 
a quantum disordered phase whose spin configuration is 
an aggregate of intra-atomic Kondo singlets. 21 ~ 25 ' When 
a magnetic field, h, is applied to the system, the mag- 
netization of this system begins to grow at the Kondo 
singlet gap. As for the Kondo necklace model without 
J2, 26 ' the magnetization curve smoothly grows with in- 
creasing the magnetic field, h, once the Kondo singlet gap 
is collapsed. 18 ' On the other hand, in the presence of the 
strong frustration (J2/J1 > 0.3419), the magnetization 
plateau appears at half saturation (m = 1/2) besides the 
Kondo singlet gap. 19 ' 20 ' The m = 1/2 plateau can be ex- 
plained as follows. The frustrated Kondo necklace model 
with m = 1/2 can be effectively regarded as spin-1/2 
zig-zag ladder which is well known as the Majumdar- 
Ghosh model 27 ' in a special case J2/J1 = 0.5, because 
the innercore spins of the system are almost ferromagnet- 
ically aligned by the magnetic field h and their degree of 
freedom is almost frozen by h. It is well known that the 
ground state of spin-1/2 zig-zag ladder is in the quantum 
disordered singlet-dimer phase. This phase corresponds 
to the m = 1/2 plateau phase in the present system. 

Since the magnetization curve of the spin-1/2 zig-zag 
ladder exhibits the cusp singularity in the vicinity of the 
saturation magnetization, we can expect that the magne- 
tization curve of the present model exhibits a similar cusp 
singularity in the high- field region 1/2 < m < 1. The ex- 
pected cusp will be found in the high-field region with 
the use of DMRG method in the next section. Moreover, 
the DMRG calculations will find an unexpected cusp just 
above the Kondo singlet gap. An appearance mechanism 
of the low-field cusp differs from that of the high-field 
cusp, because the excitations of innercore spins play an 
essential role of the magnetization process in this region. 
The main purpose of this paper is to explain the origin 
and characteristics of the low-field cusp singularity. 

This paper is organized as follows. In Sec. 2, the Hamil- 
tonian of the frustrated Kondo necklace model will be 
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Fig. 1. Schematic representation of the frustrated Kondo neck- 
lace model (1). Small and Large circles denote the conduction- 
electron spin (S c ) and the innercore spin (Sf) respectively. 

given and its magnetization curves will be calculated by 
the DMRG method. In Sec. 3, the frustrated Kondo neck- 
lace model will be mapped onto a fermionic quasiparti- 
cle system by means of the bond-operator representation 
which can explicitly describe the triplet excitation from 
the Kondo singlet sea (Sec. 3.1). A mechanism of appear- 
ance of the low-field cusp will be explained in terms of 
the double-well shape of an energy dispersion curve of 
the triplet excitation (Sec. 3.2). Moreover we will discuss 
the critical behavior of the magnetization curve near the 
cusp singularity in the low-field region < m < 1/2 
(Sec. 3.3). In Sec. 4, we will determine the phase dia- 
gram on a plane of the Kondo interaction, Jk, and the 
frustration parameter, A = J2/J1, showing whether the 
cusp singularity appears in the magnetization curve or 
not. The final section (Sec. 5) is devoted to conclusion. 

2. Hamiltonian and magnetization curves 

The Kondo necklace model is simplified version of the 
Kondo lattice model in one dimension, assuming a nar- 
row band limit for conduction electrons in the Kondo 
lattice model. In this limit, the motion of conduction elec- 
trons is prohibited, and the system is composed of two 
different kinds of localized spins: the spins of innercore 
electrons (on the / orbits) and the conduction electrons. 
Thus, the Hamiltonian of frustrated Kondo necklace is 
described as 

L L 

H = J ^ A"~ x ff C| j ■ S Ct j +v + Jk J~] S c j ■ Sfj, (1) 

i/=l,2 

which is schematically shown in Fig. 1. Here, S c ,j are 
conduction-electron spins with S c j — 1/2, and S fj are 
innercore-electron spins with Sfj — 1/2 which are hang- 
ing from the S^-spin chain. The nearest-neighbor ex- 
change coupling J and Kondo coupling Jk are assumed 
to be antiferromagnetic (J, Jk > 0). A is a frustration 
parameter which denotes competition between nearest- 
and next-nearest-neighbor interactions in S^-spins. 

In the calculation of the magnetization, we use the 
DMRG method based on the finite-size algorithm. 17 ) The 
advantage of this method is that the large system can 
be studied, as compared with the exact-diagonalization 
method. Unlike the quantum Monte Carlo simulations, 
the DMRG method is free from the negative sign prob- 
lem and can easily calculate even the ground state in the 
magnetic field. The DMRG method is an ideal tool for 
studying the magnetization curve of the one-dimensional 
spin systems. Our DMRG calculations have been per- 




Fig. 2. The magnetization curves of the frustrated Kondo neck- 
lace model with J/Jk = 1-0, In the case of A = 0.5 (filled circle), 
two cusps and a plateau appear in the magnetization curve which 
is smooth in the absence of frustration A = 0.0 (open circle). 

formed under the open boundary condition, for the sys- 
tem length L = 60 and the number of optimal states up 
to 120. 

Figure 2 shows the magnetization curves of the frus- 
trated Kondo necklace model with A = 0.0 and 0.5 for 
fixed J I Jk — 1.0. In the absence of frustration (A = 0.0), 
the magnetization curve grows smoothly with increasing 
the magnetic field h once the Kondo singlet gap Ak col- 
lapses. This means that there is no phase transition in 
the magnetic field in the range Ak < h < h$ (hs- sat- 
uration field) and the system belongs to a universality 
class of the Tomonaga-Luttinger liquid phase. 19, 20 * ) 

In the case of A = 0.5, the Fig. 2 clearly shows the 
cusp singularities in the low-field (0 < m < 1/2) and 
high-field (1/2 < m < 1) regions, which cannot be de- 
tected by the exact-diagonalization calculations. As we 
mentioned in Sec. 1, the appearance of high-field cusp 
is expected from previous results by Okunishi et a/., 15 ) 
and can be interpreted in terms of the low-energy disper- 
sion with double- well shape of single down-spin magnon 
in the up-spin sea. However, the concept of single down- 
spin magnon is useful only in the vicinity of the satura- 
tion magnetization. Therefore it cannot be used to ex- 
plain an appearance mechanism of the low-field cusp. To 
discuss the low-field cusp, we must consider appropriate 
elementary excitations in the low-field region. 

3. Mechanism of an appearance of a cusp sin- 
gularity 

3.1 Bond- operator representation 

In this section, we discuss an origin of cusp singular- 
ity in the low-field region (0 < m < 1/2). First of all, it 
should be remarked that we have to know about the low- 
energy excitation modes of the system, in order to under- 
stand the magnetization process in the low-field region. 
As was stressed before, the ground state of the present 
system (1) in the absence of the magnetic field is always 
in the Kondo singlet phase, when the Kondo coupling, 
Jk, of any strength is present. This suggests that the 



J. Phys. Soc. Jpn. 



Full Paper 



Takahiro Yamamoto, Rumi Manago, Yoshihiro Mori and Chiks 



low-energy excitations will be well described by triplet 
excitations from the Kondo singlet sea. To express this 
picture in an explicit form, we introduce the following 
bond-operators creating singlet and triplet states out of 
the vacuum |0): 



*t |0> = 


^=(IW-IW), 


(2) 


4 |o> = - 


•^(IW-IW). 


(3) 


4lo) = 


^(IW + IW), 


(4) 


A 10} = 


^(IW + IW). 


(5) 



These bond-operators were first proposed by Sachdev 
and Bhatt 28 ' in order to treat the quantum phase 
transitions between dimerized and magnetically ordered 
phases of antiferromagnets. Here, t, | and -ff, JJ. denote z- 
componcnts of S c and Sf spins, respectively. To describe 
the magnetization process of the Kondo necklace model, 
it is more convenient to use another set of bond-operators 
v) and S creating the states |flt} and |!JJ.}: 

U +|0) = ^(*t+i4)l°> = IW. ( 6 ) 

rf t io}^^ f (4-<)io} = im, (7) 

instead of the transverse triplet operators t x and t y . 29 "> 

In terms of these bond-operator representation, we can 
imagine that each lattice point of the Kondo necklace 
model is in one of the four atomic states s< |0), |0), 
t\ |0), and S |0), or equivalently, each site is occupied by 
one of the virtual particles created by bond-operators s<, 
ut, 4 an d d) . Then, these operators should satisfy the 
local constraint, 

4*i + u \ u i + *L-**J + d ) d i = *> ( 8 ) 
in order to work in the physical subspace with a single 
particle on each site. The spin operators S p j (p = c, /) 
are represented in terms of these bond-operators as 

S+j = ±= (±a)di t u] Sj + tl d dj + u]t zJ ) , (9) 
S-j = (±d] Sj T s]uj + d]t zJ + tt tjUj ) ,(10) 

S;,j = \(±s]t^±tl jSj -d]d J+ u] Ul ), (11) 

where the signs in the upper and lower rows denote the 
cases of p = c and /, respectively. 

Thus far, we have not specified the algebra of the bond- 
operators. As a matter of fact, the spin algebra: 

can be derived by using the expressions (9)-(ll), irre- 
spective that the bond-operators obey either bosonic or 
fermionic commutation relations. It is usual to assume 
for the bond-operators to be bosonic, in order to per- 
form a physically meaningful mean-field approximation 



on the basis of Bose-Einstein condensation of singlet 
states, (s) = (st) ^ 0, (t) = = 0. The bond- 

operator mean-field theory was successfully applied to 
the study of disordered dimerized phases in various quan- 
tum spin systems (e.g. spin-1/2 ladder, 30 ' Bilayer Heisen- 
berg model 31 ' and standard Kondo necklace model 21 '). 
Recently, the method of bosonic bond-operators was ap- 
plied to the studies of magnetization process of spin-1 an- 
tifcrromagnetic Heisenberg chain and spin-1/2 and spin- 
1 Heisenberg ladders, mapping onto spin-1/2 pseudospin 
Hamiltonian in the effective magnetic field. 29 ' 

In closing this section, we note that the algebra of the 
bond-operators is not important in the following discus- 
sion. An important point in this work is that the triplet 
excitations from the Kondo singlet sea are represented 
as fermionic particles, as shown in next section. Thus, 
we impose here the bosonic commutation relation on the 
bond-operators, only for convenience. Of course, even if 
we impose the fermionic commutation relation to the 
bond-operators, our results in this paper are not changed. 

3.2 Effective band theory 

We are now in the position to state a reasonable pic- 
ture of the origin of the cusp singularity in the low-field 
region, < m < 1/2, on the basis of an elementary band 
theory constructed with the use of bond-operator rep- 
resentation. To do this, the triplet excitations from the 
Kondo singlet sea are treated as fermionic particles. 

Since the ground state is always in the Kondo singlet 
phase as mentioned before, it seems to be useful to start 
from the strong Kondo-coupling limit, J/Jk <C 1. In 
the limit of strong- Jk and h ~ Ajf, we can project out 
the states t\ |0) and d^ |0), preserving only two low- lying 
states |0) and v) |0). Then, the constraint (8) can be 
rewritten as 

s]sj +u\uj = 1, (13) 
and the spin operators in Eqs. (9)-(ll) are simplified as 

Si, = T^«> jST ^7j, d4) 
S-j - T-j=s}«j = T-j=V, (15) 

S pJ = \ u ) u j = \^H- (16) 

The new operator 7T = u^Sj creates a triplet particle 
on j-th site from the Kondo singlet sea, satisfying the 
following hardcore bosonic commutation relations 

{7^7]} = 7j7j+7j7j = l, (17) 

(l]f = 7? = 0, (18) 

[li, 7]] = [71, 7j] = [li, lj\ =0, ijt j. (19) 
2 ) + 

Note that the operatots and 7^ are neither fermionic 
nor bosonic, and thanks to Eq. (17), the constraint 
(13) is satisfied automatically. The hardcore bosons can 
be mapped onto spinless fermions through the Jordan- 
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Wigner transformation, 







7 ■ exp 



exp 



' 3-1 

i?r E 7 « 7 " 

n=l 



3-1 

E 

n=l 



7j- 



(20) 



(21) 



It is convenient to introduce the Fourier transform of 
fcrmionic operator, 



(22) 



where the wave number fc takes discrete values k = 
2im/ L (n = 0, 1, • • • , L — 1) assuming L to be even in- 
teger. Using the fermionic operators and £fc, the frus- 
trated Kondo necklace model [Eq. (1)] in the magnetic 
field h is mapped onto a fermionic quasiparticle system: 

U = Eo+^(e k -h)zU k 

k 

k,k\q 

where Eq = —3JkL/4 is an energy of the Kondo singlet 
sea, 



|vac) =s|4---4l°): 



(24) 



satisfying a relation £fc |vac) = and the single-particle 
energy dispersion is given as 

e fc = J K + — [cos k + A cos (2k)} . (25) 

The last term in Eq. (23) describes the scattering be- 
tween the quasiparticles. This is negligible in the low- 
density region J2k (Xk^k) *C 1; corresponding to the 
low-magnetization m <C 1. 

We discuss the single-particle energy dispersion (25) 
in some detail. The dashed curves in Fig. 3 show the 
single-particle energy dispersion e&/ 'Jk for A = 0.2 and 
0.4 where J/Jk — 0.1 is assumed. As is seen in Fig. 3 
and Eq. (25), the shape of dispersion curve changes at the 
critical frustration A c = 1/4. (Taking account of higher 
order corrections, this value of the critical point A c = 1/4 
will be modified later.) For A < 1/4, the dispersion has 
a single bottom at k = it. On the other hand, for A > 
1/4, the dispersion curve has double- well shape with two 
bottoms at the wave numver 7r ± sin (— 1/4A) (0 < k < 
2tt) and at the energy e b /J K = 1 - {J/2J K )(X + 1/8A), 
The state at k = it, being global minimum for A < 1/4, 
turns out to be local maximum with the energy e m / Jk = 
l-(J/2Jjf)(A-l). 

These results can be also obtained by the perturbation 
expansion up to the first order in J/ Jk, and provide a 
sufficient basis to interpret the origin of the appearance 
of the cusp qualitatively. More detailed analysis up to 
the second order in J j Jk is described in the Appendix, 
and the results are shown by solid curves in Fig. 3. As 
is seen from these plots, the second order results are in 
excellent agreement with those of exact diagonalization 
indicated by open circles. 




0.95 



J/Jk=0.1 
X =0.4 



Fig. 3. The energy dispersion of a single triplet fermion for 
J/J K = 0.1 with A = 0.2 (above) and 0.4 (below). Short-dash 
curve shows the dispersion given by Eq. (25). Solid curve repre- 
sent the dispersion obtained by the perturbation expansion up to 
second order in J/Jk ( see Appendix). The open circles show the 
exact-diagonalization data at k = 2-nn/ L (n = 0, 1, • • • , L — 1) 
for L = 8, 10 and 12. 



The eigenstate with the energy efc is given as a single 
quasiparticle state, 



\k) = £ fe |vac) 



VI 



L 

E 

.7=1 



exp (ikj) s\sl 



10} • (26) 



To study the low-energy properties, one often makes 
the assumption that additional many quasiparticle eigen- 
states can be constructed as a collection of independent 
quasiparticles, that is \ kx, fe, • • • , k N ) = • • • £l N |0) 

with the energy e fcl + e k2 H efc N • 

We have prepared now to interpret the low-field mag- 
netization process just above the Kondo singlet gap 
in terms of elementary band theory with the use of above 
single-particle energy dispersion relation (25). Thereby, 
the magnetic field h plays a role of the chemical poten- 
tial or the Fermi energy and the magnetization m cor- 
responds to the quasiparticle number per site occupying 
the state below the Fermi energy h. Thus the mag- 
netization at zero temperature is given as 



1 



E 

k\ <Ckir 



i 

2^ 



d(h-e k )dk, (27) 



where kp and 6(x) are the Fermi wave-number and the 
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Fig. 4. Schematic diagram showing the relation between the double-well shaped low-energy dispersion of triplet quasiparticle and the 
cusp singularity in the low-field magnetization curve. 'G.S.' denotes the energy of ground state. 



Heaviside's step function, respectively. 

On the basis of Eq. (27), we can explain the charac- 
teristic features of magnetization curve in the low-field 
region, shown in Fig. 4. When a magnetic field h is ap- 
plied to the system, no magnetization appears until h 
arrives at the bottoms of the dispersion. If the frustra- 
tion parameter A is larger than the critical value A c in the 
region e m > h > ej, the Fermi level crosses the dispersion 
at four wave-numbers, satisfying = h, and the mag- 
netization curve shows sharp rise. With further increase 
in the magnetic field, the Fermi level eventually arrives 
at the local maximum e m at k — n, and the number 
of crossing points of the Fermi level and the dispersion 
curve suddenly reduces from 4 to 2, resulting in a sud- 
den slowdown of the rate of band filling and growth rate 
of the magnetization. This sudden change in the rate of 
band filling is reflected as the cusp singularity on the 
magnetization curve of the Kondo necklace model. 

3.3 Critical behavior near the cusp singularity 

To discuss the detailed behavior of the magnetiza- 
tion curve in the vicinity of cusp singularity, we consider 
the field derivative of the magnetization dm/dh. From 
Eq. (27), we find that dm/dh is given by the quasiparti- 
cle density of states at the Fermi energy h as follows. 



dm 
~dh 



1 

2tt 
1 

2^ 



S(h — ek)dk 



o 



de k 



dk 



(28) 



Here, J2{ tk =h} denotes a sum over all Fermi points. 
When the energy dispersion is shaped like the double- 
well, two different situations may occur (see Fig. 5(a)). 

[1] For h > e, m , there are two Fermi points ±fcp, and 
the energy dispersion can be linerized around the Fermi 
points as far as low-energy excitations are concerned: 
efc = ±vp{(k — 7r) =p kp} + h, where Dp is the Fermi 
velocity of the ^-fermions. Thus, from Eq. (28), we ob- 
tain dm/dh = L/nvp which is almost independent of h 




0.65 



0.70 0.75 

h/j K 



0.80 



Fig. 5. (a) The effective low-energy dispersion. The dashed lines 
show the magnetic field h corresponding to the Fermi level, (b) 
The field derivative of the magnetization dm/dh near the criti- 
cal point h = e m of the frustrated Kondo necklace model with 
J/Jk = 1-0 arl d ^ = 0-5- The solid line denotes asymptotic curve 
dm/dh ~ (e m — /i)" 1 / 2 near h = e m . 



because Vp only depends weakly on h. 

[2] For h < e m , there are four Fermi points satisfying 
the relation, 2/cfi — 2kp 2 = m. In this case, low-energy 
states are effectively described by the two-band model 
[See Fig. 5(a)]: 



±v F i{(k - 7r) =F k F i} 
-C(k-ir) 2 + e m , 



h, 



(29) 
(30) 
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where C is a positive constant. The parabolic dependence 
of ejp around k = ir gives decisive difference in dm/dh 
from the case of [1]. Substituting Eqs. (29) and (30) for 
Eq. (28), we obtain the most dominant contribution 

dm L 1 

dh 2nVC Ve 171 1 1 



(31) 



which diverges at h = e m . Thus, we find that the critical 
exponent characterizing the intensity of the divergence 
of the field derivative of the magnetization, dm/dh, is 
given as 1/2. 

Plots in Fig. 5(b) show dm/dh of the frustrated Kondo 
necklace model (1) with J/Jk = 1.0 and A = 0.5, ob- 
tained by numerically differentiating the DMRG data in 
Fig. 2. As is expected from the above results, for h > e rn , 
the plots depend on h very weakly. On the other hand, for 
h < e m , it diverges as dm/dh <~ (e m — h)^ 1 / 2 as h — > e m 
(the solid curve in Fig. 5(b)). Away from the cusp sin- 
gular point h — e m , the data are almost independent 
of h. This means that the low-energy dispersion is de- 
scribed by two linear-bands involving four Fermi points. 
From the above results, we find that the magnetization 
behaves as 



m - m cusp oc 



h - e r 

— V Era " 



h 



h > e r , 
h < e Ti 



(32) 



in the vicinity of the critical point h — e m with m — 
n^cusp- This result well explains the behavior of magne- 
tization curves near the cusp singularity obtained by the 
DMRG calculations (See Fig.2). 

4. Phase Diagram 

In the previous sections, we saw that the appearance 
of low-field cusp singularity in the magnetization curve 
is governed by the competition between the Kondo cou- 
pling J/Jk and frustration parameter A. In the strong 
Kondo coupling limit J/Jk —> 0, the cusp singularity 
appears when A exceeds the value A c = 1/4. In the pres- 
ence of Kondo coupling of general strength, the value of 
A c will depend on the strength of Kondo coupling. In the 
following, we determine a phase boundary on the plane 
of J/Jk and A, dividing whether the cusp appears in 
the magnetization curve or not. To do this, we note that 
the value of energy dispersion curve at zone-center point 
k = tt turns from global minimum to local maximum. 
Therefore, the critical point A c is to be determined from 
the condition that the point k — tt is the deflection point 
of the energy dispersion: 



dk 2 



= 0. 



(33) 



k — 7T 



As shown in Appendix, the critical value of frustration 
A c , including the second order correction with respect to 
J / Jk is given as 



Xr.= 



2J K + 2J 



8J K + 5 J 











32 1 


KJk) 


--( 

256 1 


\Jk) 



(34) 



In order to determine A c outside the strong Kondo- 
coupling region, however, we need the expression of en- 
ergy dispersion including the full effects of Kondo cou- 



pling and frustration. This is performed by determin- 
ing the energy dispersion for various values of J/Jk by 
means of the exact-diagonalization method. 

Here we briefly explain our exact-diagonalization 
method. Our exact-diagonalization calculations were 
performed under the periodic boundary condition S p j — 
Spj+L (p = c and /). Since the z-component of to- 
tal spin operator commutes with the frustrated Kondo 
necklace Hamiltonian (1), all eigenvalues and eigenstates 
of the Hamiltonian are labeled by the total magnetiza- 
tion M = J2j=i ( a j + Sj) ■ l n addition, since the Hamil- 
tonian (1) is invariant under the translation along the 



chain-direction (S Pt _ 



the eigenvalues and 



eigenstates are classified by the wave number k = 2irn/L 
(n = 0, 1, • • • L — 1) corresponding to eigenvalue of the 
translation operator. Thus, in Fig. 3, we plot the energy 
dispersion, E(k, M = 1) - E(k, M = 0) with k = 2nn/L, 
for three lengths L = 8, 10, and 12. 

Since the energy dispersion is an even function around 
k = 7r, wc fit the exact-diagonalization data near k = n 
using a fitting function, 



e k = AK + ^(k-n) 2 + ^(k-nf + 



(35) 



up to eighth order and estimate fitting parameters Ak, 
a, /?,••• for various A and fixed J/Jk- The parameter 
a. = (d 2 uj/dk 2 ) k=7r thus obtained for J/Jk = 0.1, 0.3 
and 0.5 are shown by circles, squares and triangles in 
Fig. 6, respectively. The next step is to interpolate a 
in order to obtain the critical point A c satisfying the 
condition (33). We assumed a liner form of interpola- 
tion function a = A(X — A c ) in analogy with that in the 
strong Kondo-coupling limit (see Appendix). Thus wc 
can obtain the critical point A c from the condition a = 
(see Fig. 6). Finally, the phase diagram on (J/Jk,X)- 
plane is constructed, performing above series of proce- 
dures for various values of J/Jk, as shown in Fig. 7. 
The short-dash curve in Fig. 7 shows the boundary given 
by Eq. (34) and is in good agreement with the exact- 
diagonalization data in the strong Kondo-coupling re- 
gion, J/Jk <C I- On the boundary, the magnetization 
curve grows rapidly at h = Ak because the low-energy 
part of the dispersion becomes quartic. To be more ex- 
plicit, the magnetization behaves as m <~ (h — Ak) 3 ^ 4 , 
as is derived from Eq. (28). 

5. Conclusion 

In this paper, the magnetization processes of frus- 
trated Kondo necklace model have been studied by us- 
ing the density matrix renormalization group (DMRG) 
method. We have paid attention to the low-field cusp 
singularity appearing just above the Kondo singlet gap. 
To reveal an appearance mechanism of the cusp, we in- 
troduced the triplet quasiparticle excited from Kondo 
singlet sea by using the bond-operator representation of 
dimerized spins, and we have discussed the mechanism 
in terms of an elementary band theory with the use of a 
single-particle energy dispersion of the triplet quasipar- 
ticle. As a result, we have made it clear that the appear- 
ance of low-field cusp has its origin in a change of shape 
of quasiparticle energy dispersion curve, reflecting effects 
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Fig. 6. The second derivative of the energy dispersion e*, at k = 
7r. The circles, squares and triangles denote (9 2 efc/9fc 2 ) fc _ 7r for 
J/JfC = 0.1, 0.3 and 0.5, respectively. Each symbol is fitted by 
a straight line. 
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Appendix: Strong Kondo coupling expansion up 
to second order in J/Jk 

In this appendix, we present the second order correc- 
tion in J I Jk to the low-energy dispersion of the frus- 
trated Kondo necklace model. 

The Kondo coupling Hamiltonian [the second term in 
Eq. (1)] will be treated as the unperturbed one and its 
ground state is given as a direct product of the Kondo 
singlets [Eq. (24)]. On the other hand, intra-chain Hamil- 
tonian [the first term in Eq. (1)] will be treated as the 
perturbation. The first order perturbation energy van- 
ishes owing to the orthogonal condition. The second- 
order correction to the energy of ground state is calcu- 
lated by the standard perturbation formula. The pertur- 
bation energy up to the second order in J/Jk is obtained 
as 



Eo/J, 



K 



1 + ^(1 + A 2 ) 




(AT) 



To obtain the low-energy dispersion, we next calculate 
the energy correction in an excited state. We suppose 
that the low-energy physics in the strong Kondo-coupling 
region is well described by the single triplet excitation 
state (26). Following the same procedure as calculating 
the ground-state energy, we can estimate the excitation 
energy. The first and the second order energies are given 
as 

I i - -^r- (cosfc + Acos(2fc)) 



Eu/Jk = 



4 



2J K 



Fig. 7. The phase diagram on plane of J/ Jk and A, which decides 
whether a cusp appears in the region of0<m<l/2in the 
magnetization curve. The short-dash curve is the boundary given 
by Eq. (34) obtained from the strong Kondo-coupling expansion 
up to second order. 



of frustration. 

We also discussed the critical behavior of the mag- 
netization in the vicinity of the cusp singularity in the 
low-field region. Calculating the density of states from 
the quasiparticle energy dispersion, we found that the 
magnetization curve behaves as m — m cusp cx — \Jt m — h 
(h < e m ) and m - m cusp cx h - e m (h > e m ). This 
result well explains the behavior of the magnetization 
curves near the cusp obtained by the DMRG calcula- 
tions. Moreover, we obtained the phase diagram on the 
( J j Jk, A)-plane deciding whether the cusp appears in 
the magnetization curve or not. 

The effective band theory based on the bond-operator 
formalism must be successfully applied to various quan- 
tum spin systems with dimer-gapped ground state. 32 -* 



1 

32 



J_ 
~Jk~ 



{(l + A 2 )(3L + 2) 



-4(2 - A) cos(fc) + 2(1 - 4A 2 ) cos(2fc) 

+4Acos(3fc) + 2A 2 cos(4fc)}. (A-2) 

Subtracting Eq. (A-2) from Eq. (AT), we obtain the 
single-particle energy dispersion up to the second order 
in J/Jk'- 



£k/JK 



(Ek — Eq)/ Jk 



1 

16 

J 



8Jk 
J 



3f) (1 + A2) 



cos k 



+ T^r{8\-(l-4X 2 )}—coa(2k) 



16Jk 

~K^' <os,:!/,i 



Jk 
( \J 



UJk 



cos(4A)3) 
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where the first order correction is consistent with 
Eq. (25). The coefficients Ak, a and (3 of the Taylor 
expansion [Eq. (35)] of the energy dispersion (A-3) are 
given by 

' 1 (l-A) 



a 



8J K 



5J 
4Jk 



4 +(3 + A) 



A 



J_ 

2 J K + 2 J 
~5J 



(A-4) 



,(A-5) 



= 



J 



8J K 

(48A 2 -41A + 5)(|A-6) 



Using the expression for a, we can determine an asymp- 
totic curve of boundary on the plane of J/Jk and A, 
dividing whether the cusp appears in the magnetization 
curve or not. Since the boundary, A c , is determined by 
the condition, a = 0, we can obtain the boundary form 
(34) as shown in Fig. 7. 
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